Astronomy & Astrophysics manuscript no. 20707 


©ESO 2013 


February 8, 2013 





On the functional form of the metallicity-giant planet correlation* 

A. Mortier 1 ' 2 , N.C. Santos 1,2 , S. Sousa 1 ' 3 , G. Israelian 3 ' 4 , M. Mayor 5 and S. Udry 5 



Centra de Astrofisica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 
e-mail: amortier@astro.up.pt 

Departamento de Ffsica e Astronomia, Faculdade de Ciencias, Universidade do Porto, Portugal 
Instituto de Astrofisica de Canarias, 38200 La Laguna, Tenerife, Spain 
Departamento de Astrofisica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain 
Observatoire de Geneve, Universite de Geneve, 51 Ch. des Maillettes, 1290 Sauverny, Switzerland 



m 
o 

(N 



CM 

OH 



Received 09/11/2012; Accepted 28/01/2013 



ABSTRACT 



Context. It is generally accepted that the presence of a giant planet is strongly dependent on the stellar metallicity. A stellar mass 
dependence has also been investigated, but this dependence does not seem as strong as the metallicity dependence. Even for metallicity, 
however, the exact form of the correlation has not been established. 

Aims. In this paper, we test several scenarios for describing the frequency of giant planets as a function of its host parameters. We 
perform this test on two volume-limited samples (from CORALIE and HARPS). 
Methods. By using a Bayesian analysis, we quantitatively compared the different scenarios. 

Results. We confirm that giant planet frequency is indeed a function of metallicity. However, there is no statistical difference between 
a constant or an exponential function for stars with subsolar metallicities contrary to what has been previously stated in the literature. 
The dependence on stellar mass could neither be confirmed nor be discarded. 

Key words. Planet-star interactions - Stars: abundances - Methods: statistical 



£h , 1. Introduction 

r^j Since the discovery of the fi rst extrasolar planet aro und a solar- 

1' like star in 1995 (51 Peg b, iMavor & Oueiozl [l995). the search 

for extrasolar planetary systems has accelerated. As of today, 
■ more than 800 planets have been announced. Most of them were 
i \ detected using the radial velocity technique. Although 800 is a 
relatively high number, th e theory of planet formation and evolu- 
tion is still being deba ted dPollack et alJl996tlMaver et al.l2 002; 
iMordasini et al.ll2009T) . The situation is particularly difficult for 
giant planet formation. C urrently, there are two proposed mod- 
els: core accretion (e.g. JPollack et al.l 1 19961 iRice & Arm itage 
120031: lAlibert et al.l [2004L where gas from the protoplanetary 
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disk is accreted around a previou sly formed rocky/icy core; an d 
the disk instability model (e.g. iBossI [19971: iMaver et alj|2002l) . 
where a planet is formed because of a direct gravitational insta- 
bility in the protoplanetary disk, in the same way as stars form 
from inte rstellar clouds. A h elpful overview of both models is 
given bv lMatsuo et al.l (120071) . 

This problem may be solved with the help of planet- 
host stars. Observational and theoretical evidence shows that 
the presence of a planet seems to depend o n several stel- 
lar pr operties, such as mass and metallicity (lUdrv & Santosl 
120071) . Concerning metallicity, it has been well-established that 
more metal-rich stars have a higher probability of harbor- 
ing a giant planet than their lower metal l icity counterparts 



( Gonzalez 11997b ISantos et all l200ll 12001: iFischer & Valenti 



20051: lUdrv & Santosl 120071: ISousa et al.l 1201 lbT 
2012). The occurrence rate even increases dramatically with in 
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* Tables 1 and 2 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



creasing metallicity. Current numbers, based on the CORALIE 
and HARPS samples, suggest that around 25% of the stars 
with twice the metal content of our Sun are orbited by a gi- 
ant plan et. This number decreases to ~ 5% for solar-metallicity 
objects (ISousa et al.ll201 lbblMavor et alJl201 ll). A similar trend 
has also been found by previous results (e.g . ISantos et al.l l2004; 
Fisc her & Valentil2 005: Johns on et al . 2010). Curiously, no such 
trend is observed for the lower mass planets (lUdrv et al. 1 12006: 
ISousa et ai1l2"0 08; May or et al.ll201 ll) . The Neptune-mass plan- 
ets found so far seem to have a rather flat me tallicity distribution 
(ISousa et al . 2008.1201 lbtlMavor et alJl20lTh . 

While this exponential trend is very clear for solar and super- 
solar metallicitie s, the situation for lo wer (subsolar) metallicities 
is still uncertain. ISantos et al. | J200l and lUdrv & Santosl (120071) 
suggest that a constant frequency may be a better fit for sub- 
solar metalli cities than a continuous e xponential. They used the 
sample from lFischer & Valentil (120051) as well as the CORALIE 
sampl e. Their suggestion has been discarded by iJohnson et al.l 
(2010), based on a Bayesian analysis on a large SPOCS sam- 
ple. Recently, it has also been suggested that there may be a 
lower limit below which no giant planets can be formed any- 
more dMortier et al.ll2.012h . The issue of giant planet frequency 
dependence on stellar properties is still not resolved. 

The observed metallicity correlation in the metal-rich re- 
gion fav ors the core-accretion model for the formation of giant 
plane ts (llda & LinlfcOOi lUdrv & Santoll2007L IMordasini et ail 
120 121) because the higher the grain content of the disk, the easier 
it is to build the cores that will later accrete gas. According to 
the disk-instability model, however, the presence of giant plan- 
ets w ould not be strongly dependent on stellar metallicity (Boss 
|2002|) . It is thus important that we fully understand what hap- 
pens around stars with different metallicities and masses since it 
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Table 1. Stellar parameters for the stars in the CORALIE sam- 
ple. The complete table is provided in electronic form only. 



Star 


[Fe/H] 


M, 


T ef f 


planet? 






(Ma) 


(K) 




HD967 


-0.68 


0.8 


5557.82 




HD1108 


0.04 


0.94 


5656.61 




HD1237 


0.07 


0.92 


5506.26 


yes 


HD1320 


-0.27 


0.89 


5688.25 




HD1388 


-0.01 


1.05 


5967.46 




HD1461 


0.19 


1.04 


5780.9 


yes 



can provide clues to the processes of planet formation and evo- 
lution. Different formation mechanisms may play a role around 
stars with different metallicities. 

In the light of this ongoing debate, we analyze two volume- 
limited samples in this paper to find out how the situation looks 
in this low-metallicity end. In Section [2] an overview is given 
of the samples and their data. Section [3] reports on the analysis 
of these samples. Finally, conclusions are made in Section |4] 
together with a discussion. 

2. The samples 

To test the giant planet frequency dependence on the stellar prop- 
erties, we used two volume-limited samples. In each sample, 
there are several planet hosts. Since this paper is about giant- 
planet frequency, we considered planet hosts as only those stars 
with at least one planet with a mass between 0.1 and 25 Jupiter 
mass. These plane ts have periods between 5 and 5000 days. 
iMavor et al.l d201 lb calculated the detection limits on these sam- 
ples. They show that the maximum orbital period out to which 
a 0.1 Jupiter-mass planet could be detected lies around 1000 - 
2000 days. 

2. 1 . Coralie sample 

The CORALIE sample is a large one of 1216 stars. It was con- 
structed from the larger CORALIE sample that consisted of 
about 1650 stars dUdrv et al.l 120001) . A color cut-off was then 
made. Stars with a B - V > 1 .2 were discarded from the sam- 
ple. This leaves us with a sampl e of about ~ 1250 stars, as used 
in seve ral previous works (e.g. ISantos et"aLll2004t [Sousa et al. 
l20TTbh . 

Since high-resolution spectra were not available for all of 
these stars, metallicities were derived using the cross correla- 
tion function (CCF ) calibration. From the Hipparcos catalog 
(Ivan Leeuwenl 12007). the color B — V was used to derive tem- 
peratu res with the calibration formula, reported in ISousa et ail 
(120081) . S tellar masses were e stimated following the same proce- 
dure as in lSousa et al.l (1201 lal) where stellar evolutionary models 
from the Padova group were applied using their web interface . 
This could not be done for 37 stars, leaving us with a final sam- 
ple of 1216 stars. The parameters can be found in Table [TJ In 
the bottom panel of Fig.Q] the data is presented in a metallicity- 
mass diagram. Metallicity and mass for the stars in this sample 
have typical error bars of 0.07 dex and 0.1 M , respectively. 

2.2. Harps sample 

Our second sample was constructed as part of a HARPS GTO 
program that aims to detect and obtain accurate orbital elements 
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Stellar mass M t 

Fig. 1. Mass - metallicity coverage of the CORALIE (top panel) 
and the HARPS (bottom panel) sample. 



of Jupiter-mass planets in a well-defined volu me of the solar 
neigh borhood (out to 57.5 pc from the Sun - lLo Curto et al.l 
2010). It can be seen as an extension to the CORALIE sample. 

Hi gh-resolution spect ra were taken with the HARPS spectro- 
graph (IMavor et al.ll2003h . In total, 582 stars had good enough 
spectra to spectroscopically derive accurate stellar parameters. 
The metalliciti es for these stars were all derived with the method 
described in lSantos et al.l (120041) . The equivalent widths were au- 
tomatically measured with the ARES code (Automat ic Routine 
for line Equivalent width s in stellar Spectra - ISousa et all 
l2007h . lScMsa et al] d201 lbh compared metallicities, derived with 
this method and the CCf correlation, for a subsample of the 
CORALIE sample. The values are consistent with a mean dif- 
ference of 0.01 dex and a dispersion of 0.07 dex. St ellar masses 
were e stimated following the same procedure as in Sous a et al.l 
(2011a), where stellar evolutionary models from the Padova 
group were applied using their web interface . All p arameters 
have been previously reported in Sousa et al. 1 (EoTTbb and can 
also be found in Table [2] In the top panel of Fig.[TJ the data is 
presented in a metallicity-mass diagram. Metallicity and mass 
for the stars in this sample have typical error bars of 0.03 dex 
and 0.025 M , respectively. 
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Table 2. Stellar parameters for the stars in the HARPS sample. 
The complete table is provided in electronic form only. 



Star 


[Fe/H] 


M, 


planet? 






(M e ) 




HD5388 


-0.28 


1.22 




HD6718 


-0.07 


0.95 


yes 


HD8038 


0.15 


1.01 




HD8535 


0.04 


1.15 


yes 


HD8930 


-0.23 


0.9 




HD8985 


-0.01 


1.2 





3. Bayesian analysis 

The analysi s performed in this pap er i s based on the methods 
described in Kass & Raftery (1995) and Johnson et a i1 d2010h .lt 
uses a Bayesian inference to find the best relationship between 
stellar mass and metallicity and the presence of a giant planet. 



3.1. Technique 

Giant-planet frequency can be expressed in terms of metal- 
licity and/or stellar mass with a particular functional form 
f(M t ,[Fe/H]), having a set of free param eters X. According 
to the theorem of Bayes (Baves & Pricd!l763l) ,we have 



P(X\d)P(d) = P(d\X)P(X) 



(1) 



where d represents the data. In our case, d can only have two 
results: a star does have a detectable planet or it does not. For 
a specific functional form /, the best set of parameters to fit the 
data will be found if P(X\d) is maximized. This means that, given 
the data, the probability that the parameters X fit the data, is at 
its maximum. 

The probability of the data itself, P(d), is constant. For the 
parameter probability, P(X), we have chosen to adopt a uniform 
distribution, allowing us to keep an open mind towards the re- 
sults. Since these two terms are constant, this means that P(X\d) 
will be maximized if P(d\X) is maximized. 

Now, assume you have N stars in your sample of which H 
stars are planet hosts. Then you can write P(d\X) as the product 
of all the separate probabilities per star: 



P(d\X) = f]P(</,|X) 

I 

H 

= Y]f(M t ,[Fe/H])[][l - f(M t ,[Fe/H])]. (2) 



N-H 



If a star has a detectable planet, its probabilty thus equals /. In 
the absence of a giant planet around the star, the probabilty is 
then obviously 1 - /. 

Since the stellar mass M* and metallicity [Fe/H] are mea- 
sured values with errorbars, it is more correct to use a proba- 
bility distribution function (pdf) p(M,-, F;) instead of their actual 
values. Per star, the giant planet frequency is then actually cal- 
culated by 



f(Mi,Fd 



II 



p{M u Fj) f(M, F)dMdF. 



(3) 



The pdf p can be approximated by a product of Gaussians with 
means (M,, F;) and standard deviations (cr M , cr fj ). 



Different functional for ms can be evaluated b y using the 
Bayesfactor, as described in lKass & Rafterv| (1 19951) . When test- 
ing two functional forms / and g, the Bayesfactor is expressed 
as 



B 



fg 



P(d\f) _ [ P(d\X f )P{Xf)dX f 
P(d\g) ~ J P(d\X g )P(X g )dXg 



(4) 



with Xf,X g the set of parameters for / and g, respectively. The 
integration limits ar e the limits of th e explored parameter space. 
According to lKass & Raftervl d!995l) . a Bayesfactor of about 100 
is needed to statistically rule out a specific functional form. 

3.2. Functional forms 



Previ o us studies (e.g. Fischer & Valenti 20051: Loyis & Mavoi 
2007; lUdrv & Santos! 12007b iJohnson et all l2010t ISousa et all 
1201 lbl) have shown that giant planet frequency is most likely ex- 
ponentially dependent on metallicity [Fe/H] and polynomially 
dependent on stellar mass M». By including those two depen- 
dencies, we adopted the following general functional form: 



f(M .\Fe/H\) = c\—\ 10 i[Fe/H] . 



(5) 



As stated before, there are still questions regarding the cor- 
rect form in the low-metallicity tail. While it must be noted that 
there are still many other functional forms that may fit the data, 
we tested seven different functional forms in our analysis: 



The traditional metallicity and mass dependence (Eq. |5}, 
Only a metallicity dependence (a = 0), 
Only a mass dependence: / ■ Ml + h, 
Eq.|5]with a fixed to 1, 

A combination of several functional forms, like 



f(M t ,[Fe/H]) = 

Eq.|6]with e = -oo, 
Eq.|6]with d — e. 



cM a JQ hlFe/H] if [Fe/H] >d 
cl0 bd if e < [Fe/H] <d 

if [Fe/H] < e 



(6) 



Form number |4] takes both metallicity and stellar mass into 
account, but fixes the mass dependence to a linear dependence. 
The last three forms will vary the behavior of the function for the 
lower metallicities. It is considered that for decreasing metallic- 
ity, the function flattens out, resulting in a constant rather than 
an exponential (paramet er d). Anoth er consideration is a sudden 
lack of giant planets. Mortie ret al.1 d2012l) suggested that there 
may be a lower limit in metallicity below which no giant planets 
can be formed. This lower limit is represented by parameter e. 
Option[5]takes both the constant and the stop into account, while 
options|6]and|7]take only the constant, resp. stop into account. 

In our parameter space, we assumed a uniform distribution 
over several intervals as seen in Table [3] These assumptions 
are based on what has been written in t he literature in previ- 
ous studies (e.g. lLovis & Mavorll2007t lUdrv & Santos! 12007b 
IJohnson et al.ll201ol) . and the derived parameters do not depend 
on the a priori chosen intervals. 

3.3. Results 

Finding the best solution of one functional form, means max- 
imizing the probability that the parameters X fit the data. 
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Table 3. Chosen interval for the parameters in the functions. 



Parameter 


Interval 


Step 


a 


[0.0, 2.0] 


0.1 


b 


[0.0,3.0] 


0.1 


c 


[0.01,0.15] 


0.01 


d 


[-0.5,0] 


0.05 


e 


[-0.75, -0.55] 


0.05 


f 


[0.05,0.55] 


0.05 


g 


[0.0,2.4] 


0.1 


h 


[-0.50, -0.05] 


0.05 



Table 4. Results of the Bayesian analysis for the three samples. 
For each sample, the functions are ordered from the best to the 
worst in the left columns. The Bayesfactors that compare two 
functions are given in the right columns. 



HARPS 


CORALIE 


HARPS+CORALIE 


Funct. 


B fg 


Funct. 


B fs 


Funct. Bf g 


4n 








K 




1.02 




1.03 


1.06 


7( 




K 




< 














1.09 




1.17 


1.20 


K 




K 




K 














1.00 




1.01 


1.01 


K 




< 




K 














1.44 




1.25 


1.19 


K 




< 




< 














1.34 




2.07 


2.71 


< 




< 




K 














73.13 


1.32 • 10 6 


> 10 10 


V 




Y 







Since there are many local maxima throughout our parameter 
space, we chose to adopt the conservative (but time-consuming) 
method of stepping through a dense grid. The step sizes we used 
are given in Table [3] To transform the products of Equation 
2] int o sums, we a lso worked in logarithmic space (see e.g. 
Johnson et al]|2010l) . 

Comparing the seven different functional forms, leads to in- 
teresting results. In Table |4] the order of functional forms is 
given for both samples, as well as the combined sample (HARPS 
+ CORALIE). The top functional form is the best, the bottom 
one the worst. On the right of these functions, the Bayesfactors 
are represented that compare these two functions. As can be 
seen, the only functional form that can statistically be ruled out 
is form number 3, where there is only a mass dependence. The 
other six forms cannot be distinguished statistically. By multi- 
plying the Bayesfactors, we calculate that the 'best' form and 
the 'worst' form of these six function relate with a Bayesfactor 
of 2.15, 3.15, and 4.12 for the HARPS, the CO RALIE, and the 
combi ned sample, respectively. According to iKass & Raftervl 
(1995), these numbers are too low to conclude anything from 
it. 

In Fig. [2] the best solutions for three functional forms are 
shown for both our samples. The functional forms that are shown 
are the traditional exponential with a linear mass dependence 
(form nr 4), the exponential with a constant (form nr 6) and the 
exponential with a constant and a drop (form nr 5). Even from 
the figures, it is clear that it is hard to distinguish between the 
different forms in the metal-poor regime. 

The conducted planet search surveys dMavor et al.l 1201 ll) 
were not complete for all period ranges. This is particularly true 




[Fe/H] 



— exponential [Fe/H], linear mass 

- - with constant 

with drop + cte 
rzzi HARPS 




[Fe/H] 



Fig. 2. Frequency of giant planets as a function of metallicity and 
mass of the HARPS (top panel), the CORALIE (middle panel) 
and the combined (bottom panel) sample. Three different func- 
tional forms are shown: a complete exponential with linear mass 
(blue curve), an exponential and a constant (green curve), and an 
exponential, a constant plus a drop (red curve). The stellar mass 
is fixed to Af« = 1 .0M o . 



for longer periods (> 1000-2000 days), where the results are in- 
complete. About a quarter of our sample have these long-period 
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planets. By excluding them, however, we obtain comparable re- 
sults from the analysis. Our conclusions are thus not sensitive to 
the inclusions of these planets. 

4. Conclusions and discussion 

We performed a statistical Bayesian analysis to look for the 
best function to describe giant planet frequency in terms of stel- 
lar properties, specifically stellar mass and metallicity. We per- 
formed this test with seven different functional forms, based on 
three volume-limited samples. One sample was observed with 
CORALIE, the other one with HARPS. The test was also per- 
formed for the combined sample. 

The test concluded that giant planet frequency is definitely 
not a function of stellar mass alone. There was no statistically 
significant result for the other functions. The Bayesfactors, used 
to compare bet ween the funct i ons, ar e all between 1.0 and 4.12. 
According to iKass & Raftery (I 1995b . this means that the differ- 
ence is not worth more than a bare mention, so it is still unclear 
what happens exactly in the metal-poor regime. 

While the dependence on metallicity is clear, an additional 
dependence on stellar mass could neither be confirmed nor dis- 
carded. The difference between the functional forms was mostly 
focused on the metal-poor regime. We tested whether the giant 
planet frequency in that regime is flat rather than exponential. 
No statistical difference was found between th e two options. 
These results contrast strongly with the results of LFohnson et alj 
d20"Tch. For our an alysis, we used the formalism as described 
by Kass & Raftervl dl995h where Bayesfactors are calculated 



as can be seen in o ur Eq. [4] and their Eqs. 1-2. Even though 
iJohnson et al.l (1201 Ot) cite the same work, it seems that they use a 
different formulation for comparing models (Eq. 7 in their work) . 
These different formulas may explain the different results. 

It has to be noted, however, that their sample has a much 
wider mass range, including both A and M stars. This makes 
their analysis more sensitive to a possible mass dependence. 
The metallicity range of their sample, however, is comparable 
to ours. 

The very clear trend that giant planet frequency is an in- 
creasing function of metallicity in the metal-rich regime can be 
explained with the theory of core-accretion. The lack of such 
a trend can be explained more easily with the theory of grav- 
itational instability. In re cent works (e.g. iMeru & Batdl2010t 
Rogers & Wadslev 2012), simulations have shown that planet 
formation induced by gravitational instability is more likely to 
occur around metal-poor stars than around metal-rich ones. 

Disk instability is expected to occur only in the outer re- 
gions, explaining among other things the e xistence of massiv e 
long-period planets like HR8799b, c and d (iMarois et al.ll2008h . 
However, this does not exclude that sho rt period plane t s may 
also be formed by gravitational instability. Barut eau et al.l d201 lb 
suggest that planets formed by gravitational instability may mi- 
grate very rapidly inwards. Several giant planets around metal- 
poor stars have periods between 50 and 1000 days (not very 
short, not very long). They may also be the result of a gravi- 
tationally unstable disk. 

These theoretical works suggest that a different trend in gi- 
ant planet frequency may be expected. If gravitational instability 
dominates in the metal-poor regime and core-accretion in the 
metal-rich regime, we would expect to see a constant followed 
by an exponential. It has to be noted that the samples that we are 
testing are simply too small to be able to distinguish between a 
constant or an exponential in the low-metallicity end. We need 
a much bigger unbiased, volume-limited sample to be able to 



form any distinctive conclusions. A sample of about 5000 stars 
will enable us to better distinguish between the different models. 
With all the current efforts to discover planets in large samples 
of solar-type stars like GAIA, Kepler, etc., it is likely that this 
problem will be solved in the near future. 
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